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ABSTRACT 

We present the resuUs of combined NH3 (1,1) and (2,2) line emission observed with the Very Large Array 
and the EfFelsberg 100 m telescope of the Infrared Dark Cloud G14. 225-0. 506. The NH3 emission reveals 
a network of filaments constituting two hub-filament systems. Hubs are associated with gas of rotational 
temperature Trot- 15 K, non-thermal velocity dispersion (Tnt~1 km s"^ and exhibit signs of star formation, 
while filaments appear to be more quiescent (Jrot-l 1 K, km s"^). Filaments are parallel in projection 

and distributed mainly along two directions, at PA~10° and 60°, and appear to be coherent in velocity. The 
averaged projected separation between adjacent filaments is between 0.5 pc and 1 pc, and the mean width of 
filaments is 0.12 pc. Cores within filaments are separated by -0.33+0.09 pc, which is consistent with the 
predicted fragmentation of an isothermal gas cylinder due to the 'sausage '-type instability. The network of 
parallel filaments observed in G14. 225-0. 506 is consistent with the gravitational instability of a thin gas layer 
threaded by magnetic fields. Overall, our data suggest that magnetic fields might play an important role in the 
alignment of filaments, and polarization measurements in the entire cloud would lend further support to this 
scenario. 

Subject headings: stars: formation — ISM: clouds — ISM: individual objects (G14. 225-0. 506) 



1. INTRODUCTION 

Filaments are ubiquitous stru ctures in star forming com- 
plexes ( e. g., [S chneider & Elme greenI IT97 9: Wiseman & Ho] 
19981: iH atchell et al. 2005; "Goldsmith et al. 2008c 
Wang et al . 2008; Jackson etal. 2010; Schneider et al. 201(1 
Molinari et all l2010l: lAndre et all 1201 Ot lArzoumanian et afl 
201 Ih . and often intersect in high-density regions of low 
aspect ratio and associated with star formation, known a s 
hub-filament systems (e.g., iMversI l2009l: iLiu et al.l l20T2h . 
However, their formation and their role in the star formation 
process is not well understood yet. 

In nearby (d -200-500 pc) molecular clouds, recent pho- 
tometric results from Herschel suggest that large-scale tur- 
bulence might be respons ible for the formation of filaments 
(lArzoumanian et al]l201 Ih . while spectroscopic studies, sen- 
sitive to smaller scales, show that filaments present subsonic 
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non-t hermal motions (iHacar & Tafallal 120111: iPineda et"an 
|2011|) . indicative of a dissipation of turbulence at smaller 
scales. While this is consistent with observations in more dis- 
tant and massive s tar-forming regions, such as G28.34-h0.06 
(IWang et al.l l2008h . a number of studies reveal supersonic 
non-thermal motions and suggest the formation of filaments 
by the converge nce of flows or by filanient-filament collisions 
on large scales (Schneider et al. "2010"; Csengeri et al." '2 0111; 
iHeit sch et al. 2008; Jimenez-Serra et al. 2010; Henshaw e t al.l 
120121: iMiettinenI 120121: iNakamuraet all 120121) . In addition, 
theoretical studies propose that magnetic fields could play a 
role in the formation of filaments (e.g., iNagai et anil998l: 
iNakamura & Lill2008h . It is clear, then, that several forma- 
tion mechanisms have been invoked to explain the formation 
and alignment of filaments and therefore further spectroscopic 
studies of filamentary regions are essential to investigate the 
origin and evolution of such structures. 

Filaments are prevailing structur es in Infrared Dark Clouds 
(IRDCs, cf. iRathbome et al.ll2006l) . In this Letter we present 
combined Very Large Array (VLA) and Eff'elsberg 100 m tele- 
scope observations of the NH3 (1,1) and (2,2) transitions to- 
ward the IRDCG14.225-0.506 (hereafter G14.2). Most of 
the studies performed so far toward this region focus on 
the brightest infrared sources, IRAS 18153-1651 (hereafter 
118153) and IRAS 18152-1658 (hereafter 118152) with a lu- 
minosity of -1.1x10"^ Lq and -4x 10^ Lq, respectiv e ly, an d 
located at a distance of 2.3 kpc (iJaff'e et al.l Il98ll [l982h . 
Single-dish observation s show that 118153 is associated with 
H2O maser emission (Ja fl^e et al.ll9 81: Pala d et al.ll993h. and 
dense gas emission (iPlume et al.l [1992: Anglada et al.lll996l: 
iBronfman et al.l [T99a) ^ More recent VLA obse rvations re- 
veal H 2O maser emission in 9 difl'erent positions (IWang et al.l 
l2006b . which indicates that star formation is already ongoing 
in some parts of the cloud. IRDC G14.2 has been identified, 
using Spitzer data, by lPeretto & Fulled (l2009h as a cloud con- 
taining an important amount (-100) of density enhancements 
or fragments displaying a filamentary morphology. The fila- 
mentary appearance of G 14. 2 (Fig.[TJ and its relatively nearby 
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Fig. 1. — Archival Spitzer 4.5/8.0/24 /um (blue/green/red) three-color com- 
posite image of IRDC G14.225-0.506. 



distance makes this region a good selection to investigate the 
physical properties of filaments and their formation mecha- 
nism at high spatial resolution. 

2. OBSERVATIONS 

The observations of the NH3 (1,1) and (2,2) transitions were 
conducted using the VLW^ in the D configuration on Novem- 
ber 12 2005 (project AW666). We performed a 34-pointing 
mosaic covering an area of 7' x 13'. The integration time was 
-4.5 minutes per pointing. The adopted flux density of the 
flux calibrator 3C 286 was 2.41 Jy at a wavelength of 1.3 cm. 
The time variation of the gains was calibrated using J1832- 
105, with a bootstrapped flux of 0.97+0.01 Jy, and the band- 
pass calibrator used was 3C 273. We used the 4IF spectral line 
mode, which allows simultaneous observations of the NH3 (1, 
1) and (2, 2) lines with two polarizations for each line. The 
bandwidth used was 3.12 MHz, divided into 63 channels with 
a channel spacing of 48.8 kHz (-0.6 km s"^ at 1.3 cm), cen- 
tered at -21 km s"^ The visibility data sets were calibrated 
using the AIPS software package of the NRAO. 

To recover extended structures filtered out by the interfer- 
ometer, we performed NH3 observations with the Eff'elsberg 
100 m telescope (project 101-07). The observations were 
carried out between April 4 and 7 2008. We used the 18- 
26 GHz HEMT receiver tuned to a frequency of 23.7 GHz 
with the 16384 channel Fast Fourier Transform Spectrome- 
ter, allowing simultaneous observations of the NH3 (1,1) and 

The Very Large Array (VLA) is operated by the National Radio As- 
tronomy Observatory (NRAO), a facility of the National Science Foundation 
operated under cooperative agreement by Associated Universities, Inc. 



(2,2) lines. The total bandwidth used was 100 MHz, which 
provides a velocity resolution of 0.075 km s"^. The obser- 
vations were conducted in frequency switching mode with a 
frequency throw of 7.5 MHz. At the observed wavelength, 
the half-power beam width of the telescope is ~40''. The 
map covered an area of 8' x 13' and was made by observ- 
ing the positions of a grid with half-beam spacing. Pointing 
was checked at hourly intervals, with a pointing accuracy bet- 
ter than 8''. To convert the arbitrary noise tube units of the 
Effelsberg data to main beam brightness temperature we ob- 
served as a primary flux calibrator NGC 7027 and a nearby 
quasar as a secondary flux calibrator. Data reduction was 
performed using the CLASS package, which is part of the 
GILDAS0 sofware. We combined the visibility data from 
the VLA and Effelsberg 100 m telescope for both NH3 (1,1) 
and N H3 (2, 2) lines follo wing the MIRIAD procedure out- 
lined in I Vo gel et al.l (119841) . We applied a wf;-taper function of 
23 k/l during imaging. The resulting synthesized beams were 
8r2x7r0 (PA= -15°) for NH3 (1,1) and 870x679 (PA= -16°) 
for NH3 (2,2). The rms was ~8 mJy beam"^ per 0.6 km s"^ 
spectral channel. 

3. RESULTS AND ANALYSIS 

Figure [2] (top left) shows the combined (VLA-hEfl'elsberg) 
zero-order moment map of the NH3 (1,1) emission overlaid 
on the 8 jim Spitzer image. The overall morphology of the 
NH3(1,1) dense gas consists of extended and clumpy fila- 
mentary structures, strikingly mimicking the extinction fea- 
ture seen in the Spitzer image. While the NH3 (1,1) emis- 
sion is spatially extended, the NH3 (2,2) emission is compact 
(Fig. [21- top right), suggesting that the extended emission is 
at lower temperatures. We identified the most prominent fil- 
aments based on the morphology of the NH3 (1,1) together 
with the fact that these structures are coherent in velocity. We 
used the following criteria: i) filaments must have aspect ratio 
larger than 6; ii) the signal-to-noise ratio should be larger than 
$3 and Hi) they must appear in at least two velocity chan- 
nels and spanning a maximum velocity range of 3 km s"^ 
Fig. [21 (top right) shows, for comparison, the 870 jim con- 
tinuum en iission from th e LAB OCA bolometer at the APEX 
telescope (lBusquetll2010l) . supporting our identification. 

We identified a network of 8 filaments and two hubs (named 
hub-N and hub-S in Fig. [21), which were recognized using 
the NH3 (2,2) emission as denser regions in which some fil- 
aments intersect. The NH3 filaments, which cover a total area 
of 4.7x8.7 pc, appear approximately parallel, in projection, in 
two preferred directions, at PA of 10° and 60°, and they con- 
tain chains of dense core£3 aligned along the filament axis 
and distributed at somewhat regular spacings of about -30'' 
or 0.33+0.09 pc at the distance of the cloud. The averaged 
projected separation between adjacent filaments is between 
0.5 pc and 1 pc. In Table [T] we report on the length and width 
at Full Width Half Maximum (FWHM) of each filament ob- 
tained from NH3 (1,1) data. On average, we found that the as- 
pect ratio is -15:1, with a typical FWHM width of -0.12 pc. 
This value is close to the filament width of 0.1 pc reported 
for the IC5146, Aquila, and Polaris molecular clouds from 

13 http://www.iram.fr/IRAMFR/GILDAS 

1^ Signal-to-noise ratio computed in the zero-order moment map of 
NH3 (1,1), where the rms noise level has been estimated using 
where cr is the rms noise of the channel maps and Av=0.6 km s"^ . 

1^ Cores have been identified in the zero-order moment map of NH3 (1,1) 
if at least the 6cr level is closed, where cr is the rms noise of the map. 
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Fig. 2. — Top Left: Combined NH3 (1,1) integrated intensity (contours) overlaid on the 8 //m Spitzer image (color scale). Contour levels range from 3 to 18 
in steps of 3, and from 18 to 58 in steps of 10 times the rms noise of the map, 9 mJy beam~^ km s~^. The dashed line indicates 50 % of the sensitivity level 
of the VLA mosaic. Gray lines are the polarization vectors of the near-infrared (//-band) observations (Busquet et al. in prep.) with the scale shown in the top 
right comer. Top Right: Combined N H3 (2,2) integra ted intensity (gray scale and black contours) overlaid on the 870 jim continuum from LAB OCA bolometer 
at the APEX telescope (gray contours; IBusquetll201QI) . Contour levels for NH3 (2,2) range from 2 to 10 in steps of 2, and from 10 to 60 in steps of 5 times the 
rms noise of the map, 9 mJy beam~^ km s~^ Contour levels for 870 jim emission range from 3 to 53 in steps of 10, and from 53 to 653 in steps of 100 times 
the rms of the map, 25 mJy beam~^. The NH3 and 870 //m continuum synthesized beams are shown in the bottom left and bottom right comers, respectively. 
Bottom: Contours: zero-order moment map of NH3 (1,1). Color scale: NH3 (1,1) first-order moment map {Left) and second-order moment map (Right). Units 
are km s~^ . Stars indicate IRAS sources in the field, and crosses mark the position of H2O maser (Wang et al. 2006). The most prominent filaments are indicated 
with white/black lines and labelled according to its position angle. Arrows in the bottom right panel mark the positions of the NH3 (1,1) spectra shown in Fig. [3] 
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Fig. 3. — NH3 (1,1) spectra, in units of brightness temperature, at some se- 
lected positions depicted with arrows in FigHU 



Herschel observations (lArzoumanian et al]l2Ql Ih . 

In Fig. [21 (bottom left) we present the first-order moment 
map of the NH3 (1,1) main line. Within each filament the ve- 
locity variations are small, in the range of 1-2 k m s~^ (see Ta- 
ble[T] ), similar to other filamentary IRDCs (e. g.. I Jackson et al.l 
l201Qh . This network of filaments seems to be separated into 
two main velocity components, one at f;LSR~ 18.3-20.8 km s"^ 
and another one at f;LSR~20.8-23.2 km s"^ which overlap in 
the hubs. The second-order moment map is presented in Fig. |2] 
(bottom right), and shows that the velocity dispersion is lo- 
cally enhanced (cr-l km s"^) toward hubs. Additionally, a 
high velocity dispersion ((T~1.6 km s"^) is seen toward an arc- 
shaped structure connecting filament FIO-E with the southern 
filaments, in a small region intersecting filament F60-C1 and 
labelled as position 'e' in Fig.[2l In this region the large values 
of the velocity dispersion are due to the presence of two ve- 
locity components separated by ~3 km s"^ (see Fig.[3]-e). The 
presence of two velocity components is also found in regions 
where filaments intersect hubs (see Fig. [3]-d,f). In contrast, 
all the other filaments appear more quiescent, with a typical 
velocity dispersion of -0.4-0.6 km s"^. 

To obtain the main physical properties (rotational temper- 
ature Trot, total velocity dispersion (Ttot, and mass per unit 
length Ml I) of each filament, we extracted an averaged spec- 
trum of NH3 (1,1) and (2,2) over the filament area at FWHM. 
The results are reported in Table [T] The rotational tempera- 
ture ranges between 10 K and 16 K. The total velocity disper- 
sion of the gas ranges from ~0.5 km s"^ up to 1.1 km s"^ and 
the non-thermal velocity dispersion over the isothermal sound 
speed, ct^t/Cs, ranges between 2 and 5, implying that fila- 
ments in G 14. 2 are characterized by supersonic non-thermal 
motions. In Table [T] we also list the mass and virial mass per 
unit length, the observed separation between cores, and the 
number of cores in each filament. Finally, the total surface 
density estimated by taking the spectrum averaged over all 
NH3 filaments is S^l.9 gcm"^. 

4. DISCUSSION AND CONCLUSIONS 

In the previous section we presented the main properties of 
the two hub-filament systems in G14. 225-0. 506. Hubs are 
more compact (aspect ratio 5 vs 15), warmer (ri-ot-15 K vs 
1 1 K), and show larger velocity dispersion and larger masses 
per unit length than filaments. Interestingly, hubs are associ- 
ated with H2O maser emission (Wang et al. 2006) and mid- 
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Fig. 4.— Top: NH3 (2,2)/NH3 (1,1) map around 118153 (color scale) over- 
laid on the NH3(1,1) integrated intensity (contours). The long red arrow 
depicts the position-velocity (PV) cut, where the center of the cut is indicated 
by the intersection of the two red lines. Bottom: PV plot of the NH3 (1,1) 
emission along the cut at PA=10°. Contours start at 0.02 and increase in 
steps of 0.02 Jy beam~^ . Positive offsets increase as indicated by the arrow. 



infrared sources (see Figs.[T]|2}, and they are the main sites of 
stellar activity within the cloud. 

The stability of the filaments ca n be studied by estimatin g 
the virial parameter ayir=Myir/M (iBertoldi & McKed[T992h . 
which is <2 for all the filaments and hubs except for FIO- 
E, and five out of eight filaments are near virial equilibrium 
(a^vir-l). This indicates that most of the filaments are unstable 
(collapsing) and probably undergoing fragmentation, compat- 
ible with the clumpy structure of G 14. 2. It is worth noting that 
filament FIO-E has T^ot and velocity dispersion values similar 
to the hub properties. Filament FIO-E presents some striations 
converging toward it. However, while hub-N and hub-S seem 
to be places where two diff'erent velocity structures converge, 
FIO-E shows only one velocity component. We speculate that 
118153, an UC HII region with L-IO"^ L© (Sanchez-Monge, 
private communication), may compress the gas, heating and 
injecting turbulence to this filament (a^vir -5). The interaction 
of this UC HII region with the dense gas is also seen in hub- 
N, where the NH3 (2,2)/NH3 (1,1) map shows a local heating 
(Fig.|4]). The position- velocity (PV) plot along this hub (see 
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TABLE 1 

Physical properties of filaments and hubs 





length^ 


width^ 


aspect 


l^LSR 




(T tot'' 


Cnt/Cs'^ 


Mil'' ^ 






/iobs^ 




A^cores^ 


Filament 


(pc) 


(pc) 


ratio 


(km s-i) 


(K) 


(km s-i) 




(Mopc-i) 


(Mopc-i) 




(pc) 


(pc) 




FIO-W 


1.48 


0.11 


13.5 


20.8-22.0 


10±2 


0.49 


2.38 


105 


112 


1.1 


0.27±0.08 


0.16-0.40 


6 


FIO-C 


1.18 


0.14 


8.4 


20.8-22.0 


10±1 


0.50 


2.31 


147 


115 


0.8 


0.33±0.07 


0.16-0.41 


4 


FIO-E 


1.63 


0.15 


10.9 


20.2-22.0 


16±1 


0.98 


3.63 


87 


448 


5.1 


0.32±0.11 


0.21-0.80 


10 


F60-N 


1.23 


0.09 


13.7 


19.5-20.8 


12±2 


0.54 


2.31 


149 


135 


0.9 


0.25±0.09 


0.18-0.44 


7 


F60-C1 


1.81 


0.12 


15.1 


18.3-20.2 


10±3 


0.63 


3.10 


120 


185 


1.5 


0.38±0.11 


0.16-0.52 


8 


F60-C2 


0.66 


0.09 


7.3 


20.8-22.6 


10±3 


0.61 


3.03 


74 


174 


2.3 


0.49±0.09 


0.16-0.50 


5 


F60-C3 


3.34 


0.12 


27.8 


21.4-23.2 


10±1 


0.73 


3.67 


218 


248 


1.1 


0.25±0.12 


0.16-0.60 


14 


F60-S 


2.12 


0.10 


21.2 


18.3-20.8 


10±3 


0.98 


5.02 


132 


447 


3.4 


0.38±0.11 


0.16-0.80 


6 


Hub-N 


1.12 


0.23 


4.9 


17.7-20.8 


15±1 


0.98 


3.80 


266 


446 


1.7 








Hub-S 


1.15 


0.23 


5.0 


18.3-20.8 


15±1 


1.09 


4.31 


328 


548 


1.7 









^ Deconvolved size at Full Width Half Maximum (FWHM) not corrected for projection effects. 



^ Trot has been derived following the appendix of lBusquet et al.l 120091) . 

^ Averaged values within the area at FWHM. 

Mass per unit length, where the mass, M = N(li2) l.Srrip A, has been estimated assuming an NH3 abundance of 3 x 10"^ (average value measured in IRDCs; Pillai et al. 2006) 
and using the area A of the filament at FWHM. The uncertainty in the mass is a factor of 3. 

^ Virial mass per unit length M^iJl - Icr^Q^jG, and virial paramter avh- = My[JM jBertoldi & McKed 19921) 

^ /lobs: observed separation between cores within a filament. Af-. predicted core separation Af = 22 H, where H = Csi^nGl.Smu «c)~^^^ is the scale height, with Cs the isothermal 
sound speed (estimated by converting Trot to kinetic temperature using the expression of lTafalla et al.l i2004f) ), G the gravitational constant, and ric the gas density at the center 
of the filament, adopted to be 10^ cm~^. The first value corresponds to the core separation using and the second value was obtained replacing by cttot- 

^ Number of cores within each filament. 



Fig. [4]) reveals an inverted C-Hke structure, consistent with an 
expanding shell (lArce et al.ll2011h . 

We investigated the fragmentation of filaments in the mag- 
netohydrodynamic 'sausage '-type instability scenario (Chan- 
drasekhar & Fermi 1953, see also Jackson et al. 2010), which 
predicts periodic separation between fragments (or cores) for 
a given density and isothermal sound speed. For an isother- 
mal gas cylinder of finite radius R, the core separation can be 
expressed as Af=22H for H, where H is the scale height 
(see Table [T] for the formal expression). This is the case of 
G14.2, since R and H are 0.12 pc and 0.04 pc, respectively. 
Adopting a density of 10^ cm"^, we estimated the predicted 
core separation using first the isothermal sound speed, yield- 
ing /If ~0. 16-0.21 pc, and then replacing Cs by the total veloc- 
ity dispersion (Ttot, which gives /If -0.4-0.8 pc (see Table [TJ. 
The observed separation, -0.33 pc, is in agreement with these 
two extreme cases. It is noteworthy that most of the cores ap- 
pear to be elongated along the direction of the filament, which 
could imply the possibility of further fragmentation at smaller 
scales as observed in the IRDC G28.34-h0.06-Pl (Zhan g et al.l 
l2009l:IWang etani20TTh . 

One of the most intriguing features of G 14. 2 is the network 
of filaments that are aligned in parallel. The filaments ap- 
pear to take two preferred directions, one group at a PA of 
10°, and the others at a PA of 60°. This network of fila- 
ments may arise from a layer of self-gravitating gas. Instabil- 
ity analysis has been performed for a layer of an isothermal 
infinite sheet under dynamical perturbation (|Ledoux|[r951l; 
|Schmid-Bur2k"1967'; ^Elmegreen & Elmegreenlll978l; iLarsonl 
[1985; Nagaietal. 1998; Myers 2009|). The gas is unstable 
to perturbations, which leads to high density columns in the 
plane as a result of gravitational instability. The spacing s 
between the high density columns of gas correspond to the 
wavelength of the fastest growth mode. In the absence of 
magnetic fields, the growth of instability does not have a pre- 
ferred direction in the plane. As a result, a grid of connected 
filaments appears in the gas layer. If the gas layer is threaded 
by magnetic fields, the growth of the instability develops un- 



restricted in one di rection and i s supp ressed along the or- 
thogonal direction. iNagai et aH (119981) analyzed a pressure 
confined isothermal gas layer threaded by uniform magnetic 
fields. They found that in the regime of smaller external pres- 
sure (i. e., the scale height if <^Zb, where Zb is the thickness of 
the gas layer) the instability grows faster along the field lines. 
As a consequence, high density columns, or filaments develop 
with their longitudinal axis perpendicular to the field lines. In 
the high pressure regime (i. e., /f»Zb), the fastest growth of 
instability is perpendicular to the field lines and gives rise to 
filaments parallel to the magnetic fields. 

Recent numerical simulations (Van L oo, private commu - 
nication) confirm the linear analysis in iNagai et aH (119981) . 
The simulations show that an array of high density columns 
develop in the gas layer with magnetic fields. In addition, 
lower density filamentary structures are also present, inter- 
connecting the main filaments. The highest density structures 
are found at the intersections of major and minor filaments, 
as in this work. Furthermore, the simulations show that grav- 
itational instability develops within a filament during the fila- 
ment formation. This means that fragmentation of a filament 
into cores occurs simultan eously with the fra gmentation of 
the sheet, but according to iToala et all (1201 2|) with diff'erent 
free-fall times. 

The array of filamentary structures in G 14. 225-0. 5 06 may 
arise from gravitational instability of a thin gas layer with 
magnetic fields. In fact, preliminary near-infrared polari met- 
ric observations around hub-N (Busquet et al. in prep.) re- 
veal that the magnetic field is perpendicular to filaments at 
PA- 60 ° (see Fig. [2]). Therefore, according to iNagai et aP 
(Il998h G14.2 would be in the regime of small external pres- 
sure (H ^ Zb). Using the total surface density (see Sect. 3), 
the scale height H of the initial gas layer is H -0.09 pc. This 
value should be regarded with caution, and to definitively as- 
sess its validity one would need observations of a low-density 
gas tracer to be sensitive to the gas layer. The wavelength of 
the fastest mode c an be expressed as /Ifastest = 4;: // (Eq. 60 in 
iNagai et anil998h . Using our estimation of //, the predicted 
separation is - 1 . 1 pc, in agreement with the observed filament 
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separation (between 0.5 and 1 pc). It is not clear how such a 
large gas layer (4.7x8.7 pc) may form initially. The conver- 
gence of dynamic flows could be responsible for the formation 
of such a large gas layer that subsequently could fragment into 
parallel filaments as a result of magnetic modulation. Our 
NH3 data, although showing two velocity components, do 
not reveal evidence of converging/interacting flows and fur- 
ther obs ervations of low-velocity shock tracers, like Si O or 
CH3CN (iJimenez-Serra et al.ll2Ql()HCsengeri et aDl2Ql Ih . and 
a tracer of low density material are required to definitely iden- 
tify signatures of converging flows. Overall, our data suggest 
that magnetic fields might play an important role in the align- 
ment of filaments, and polarization measurements in the entire 



cloud would lend further support to this scenario. 
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